Available  online  at  http://scik.org 
J.  Math.  Comput.  Sci.  3  (2013),  No.  6,  1353-1372 
ISSN:  1927-5307 

PARAMETER  ESTIMATION  FOR  A  THIN  LAYER  BY  MEASURING 
TEMPERATURE  INDUCED  BY  A  HEAT  SOURCE 

HONGYUN  WANG1,  HONG  ZHOU2’* 

’Department  of  Applied  Mathematics  and  Statistics,  University  of  California,  Santa  Cruz,  CA  95060 

2Department  of  Applied  Mathematics,  Naval  Postgraduate  School,  Monterey,  CA  93943 

Copyright  ©2013  Wang  and  Zhou.  This  is  an  open  access  article  distributed  under  the  Creative  Commons  Attribution  License,  which 
permits  unrestricted  use,  distribution,  and  reproduction  in  any  medium,  provided  the  original  work  is  properly  cited. 

Abstract.  We  study  the  temperature  profile  for  a  semi-infinite  thin  layer  due  to  a  Gaussian  heat  source.  We 
follow  the  time-dependent  temperature  profile  at  the  center  of  the  top  surface  of  the  layer,  deriving  the  asymptotic 
behavior  and  both  the  short  time  and  long  time  behaviors.  We  also  show  that  one  can  estimate  thermal  properites 
from  various  time  measurements  of  the  temperature  at  the  center  of  the  top  surface.  Finally,  we  investigate  the 
statistical  uncertainty  in  estimating  parameter  values  from  temperature  measurements. 
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1.  Introduction 

Laser  heating  is  an  important  technology  in  materials  processing  [14].  Since  1970s  there  are 
many  studies  on  the  temperature  profiles  induced  by  laser  radiation.  A  few  examples  include 
the  spatial  distribution  of  the  temperature  rise  due  to  a  stationary  Gaussain  laser  beam  in  a  solid 
[9,  10,  1 1],  time-dependent  solution  for  a  scanning  Gaussin  beam  [7,  15],  laser-melted  front  [4], 
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continuous  wave  (cw)  laser  annealing  of  heterogeneous  multilayer  structure  [6],  temperature 
profiles  due  to  a  moving  ellptical  cw  laser  beam  [12],  the  solution  of  the  transient  heat  equation 
with  mixed  boundary  conditions  [1],  the  Green’s  function  solution  for  a  two-layer  structure 
with  scanning  circular  beams  [5],  solution  for  moving  semi-infinite  medium  under  the  effect 
of  a  moving  laser  heat  source  [2],  time-dependent  heat  conduction  in  a  thin  metal  film  [17], 
solution  of  the  temperature  distribution  in  a  finite  solid  caused  by  a  moving  heat  source  [3], 
analysis  of  short-pulse  laser  heating  of  metals  [13],  temperature  rise  induced  by  a  rotating  or 
dithering  laser  beam  in  various  solid  structures  [16,  18,  19,  20], 

Despite  all  these  advances  in  the  studies  of  laser  heating,  it  is  still  desirable  to  obtain  analyti¬ 
cal  results  even  for  simple  cases.  In  this  paper  we  consider  the  temperature  profile  of  an  infinite 
thin  layer  produced  by  a  heat  source  in  the  form  of  a  stationary  Gaussian  beam  and  reveal  how  to 
estimate  thermal  properties  from  the  temperature  measurements.  We  first  give  the  mathematical 
formulation  and  then  obtain  analytical  solution,  paying  special  attention  to  the  time-dependent 
temperature  at  the  center  of  the  heat  source  at  the  top  surface.  We  further  discuss  the  asymptotic 
behaviors  of  the  temperature  at  the  center  of  the  top  surface  for  several  cases.  In  addition,  we 
show  how  to  derive  thermal  properties  from  the  measurements  of  the  temperature  profile  at  the 
center  of  the  top  surface,  including  short  time  data,  long  time  data  or  a  full  range  of  time  data. 
We  then  evaluate  the  statistical  error  in  the  parameter  estimations.  Finally,  we  summarize  our 
results. 

2.  Mathematical  formulation  and  analytical  solution 

Consider  a  thin  layer,  which  is  infinite  in  both  the  x-  and  y-  dimensions,  with  a  thickness  £  in 
the  z-dimcnsion.  We  establish  the  coordinate  system  such  that  the  top  surface  of  the  layer  is  at 
z  =  0  and  the  bottom  surface  is  at  z  —  —£■ 

We  assume  that  in  the  z-  direction,  the  heat  source  is  concentrated  at  a  level  immediately 
below  the  top  surface  (z,  =  0),  as  illustrated  in  Figure  1. 

In  the  x-y  plane,  the  heat  source  has  a  Gaussian  distribution  centered  at  the  origin  with  s- 
tandard  deviation  a  (we  shall  call  a  the  radius  of  heat  source).  Thus,  the  heat  source  has  the 
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heat  source 


Figure  1.  A  heat  source  is  concentrated  at  a  level  immediately  below  the  top 
surface  of  a  thin  layer  where  z  —  0. 

distribution  defined  by 

(1)  <t>(x,y,Z )  =  total  ■  /(*,  G2 )./(>’,  o2)8(z  -  0_) 


where  8  is  the  Dirac  Delta  function  and 

(2)  f(x'°2)=v^txp&) 

is  the  probability  density  function  of  a  Gaussian  distribution  and  (j) total  is  the  total  heat  flux  of 
the  heat  source  (total  heat  input  per  unit  time)  given  by 

(3)  had  =  f  (j)(x,y,z)dxdydz. 


Let  us  consider  the  situation  where  the  thin  layer  is  insulated  at  both  the  top  surface  and  the 
bottom  surface.  Let  u(x,y,z,t )  denote  the  temperature  at  position  (jc,y,z)  at  time  I.  Then  u  is 


governed  by  the  heat  equation  with  boundary  and  initial  conditions: 


(4) 


du 

~di 


K  f  d2u  d2u  d2u\  1 

^Xdx,-  +  3?  +  3?)+plrrHX’y’Z) 


du 

dz 


z=0 


0, 


du 

dz 


=  0 


Z=~£ 


u\t=0  ~  W(E 
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where 


•  K  is  the  thermal  conductivity, 

•  p  is  the  density,  and 

•  cp  is  the  specific  heat. 

The  boundary  conditions  of  zero  heat  flux  at  z  —  0  and  at  z  —  —£  can  be  enforced  by  extending 
function  u  in  the  z-  direction.  We  first  extend  u  to  z  G  [— e,  e]  by  mirror  reflection  with  respect 
to  z  =  0. 

u(z)  =  u(—z)  for  z  e  [0,e]. 


Then  we  extend  u  to  z  G  (— °°,  °°)  by  periodical  extension  with  period  =  2e: 


u(z  +  2e)  =  u(z). 


After  the  extension,  u  satisfies  the  initial  value  problem  (IVP): 
K 


(5) 


du 

~di 


d2u  d2 


+ 


u  d2i 

+ 


pcp  \dx2  dy2  dz2 


-|-oo 


+  £  S(z-2 je) 


Pcp 


)=- 


U\t=0  =  MQ- 


To  solve  the  IVP  (5),  we  first  carry  out  non-dimensionalization.  For  mathematical  conve¬ 
nience,  we  shall  denote  the  physical  variable  before  non-dimensionalization  by  xp\iy  and  denote 
the  variable  after  non-dimensionalization  simply  as  x  without  a  subscript.  We  nondimensional- 
ize  the  physical  variables  as  follows. 

1  1  1 

x  xphy  —  ■>  y  y  phy  ~i  Z  Zphy  ~  ■> 


f  —  ^  phy 


K 


pCp(J 


2’ 


(6) 


{ttphy  Mq)  ; 


KG 


After  non-dimensionalization,  the  initial  value  problem  (5)  becomes 


(7) 


du 

Tt 


d2u  d2u  d2u 


-|-oo 


^c2+^2+^2j+/(X’1)/^1).E  5^~2j 


0. 


m|/=o  = 
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The  fundamental  solution  of  the  one-dimensional  heat  equation  (also  called  the  heat  kernel) 
can  be  expressed  in  terms  of  the  Gaussian  density  as  [8] 

<8) 

The  solution  of  (7)  is  then  given  by 

u(x,y,z,t)  =  [  [  f{x-x!,2{t-t'))f(y-y',2{t-t'))f(z-z!,2{t-t')) 

JO  J{x',y',z') 

+°°  e 

(9)  x/(x',l)/(y',l)  Y  8(z' -2j—)dxdy'dzdt'. 

j=—°° 

Recall  the  fact  that  the  convolution  of  two  Gaussian  densities  is  still  a  Gaussian  density: 


x <*\ i  )/(-*',  of)dx  =  f{x,  of  +  of). 

Using  this  fact  and  making  a  change  of  variables  x  —  t  —  t1 ,  we  can  rewrite  the  solution  as 

ft  +°°  /  £  \ 

(10)  u(x,y,z,t)  =  /  f{x,2x  +  l)/(y,2r+  1)  Y  f  (z-2j-,2x)  dx. 

JO  jtJ  oc  v  > 


From  an  experimental  point  of  view,  we  are  interested  in  the  temperature  profile  at  the  ori¬ 
gin  (jc,y,z)  =  (0,0,0),  which  is  exactly  the  center  of  the  heat  source  at  the  top  surface.  The 
temperature  at  the  origin  has  the  expression: 

rt  1 


(ID 


w(0,0,0,  t)  = 


1  h(*f\dT, 


Jo  2k(2x  + 1)  y/Anx  \£ 


where  function  h{r\)  is  defined  as 

-|-oo 

(12)  Kt])=  Y  exP(-//T7)- 


j=- 


We  apply  another  change  of  variables  s  —  \fx  to  obtain  the  analytical  form  of  the  temperature 
at  the  center  of  the  top  surface: 


(13) 


1  ['ft  1 

w(0,0,0,t)  = — r  /  — —h 


2x2  Jo  (2v2  +  1)  V  £ 


<xV 


ds. 


Now  we  study  the  properties  of  function  h{r\)  defined  in  (12). 
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For  small  r\,h(r\)  is  very  close  to  1  since 


h(r\)  =  1  +  2exp(— 1/17)  +  2exp(— 4/17)  H - 

=  1  +  T.S.T.  (Transcendentally  Small  Term) 


For  large  r\ ,  let  Av  =  1  /  y/rf,  we  have 


-|-oo 

HV)  =  £  exp(-//T7) 

j=-oo 

-\-oo 

=  Vn  £  exp(-(y’Av)2)Av 
j=- 00 

r+°°  9 

~  y/ff  /  exp(—  x/dx 

J  — 00 


=  VvVn. 


Actually,  for  large  rj ,  we  have 


h{r])  —  y/rf  T.S.T.  (Transcendentally  Small  Term) 


The  plot  of  h(r\),  together  with  1  and  y/ffn,  is  shown  in  Figure  2.  Figure  2  confirms  the 
above  asymptotic  results.  Namely,  for  small  values  of  77,  h(rj)  can  be  closely  approximated  by 
the  constant  function  1;  for  large  values  of  77,  h(rj)  can  be  well  approximated  by  /r\K. 

3.  Asymptotic  solutions  at  the  center  of  the  top  surface 

Having  obtained  the  analytical  form  of  the  time-dependent  temperature  at  the  center  of  the 
top  surface,  namely,  m(0,0,0,?),  we  now  turn  our  attention  to  its  asymptotic  behaviors. 
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Figure  2.  The  plot  of  functions  h{r\),  1  and  y/rjn. 


For  small  t  (when  to2 /e1  «  1),  the  temperature  at  the  origin  can  be  written  as 


w(0,0,0,r)  = 


1 

rVt 

2  712 

Jo 

1 

r\ft 

2  TV- 

Jo 

1 

1 

2^v/2“ 

1 

—  _  3 

Vt( 

-h[^\ds 


-ds  +  T.S.T 


arctan(v/2t)  +  T.S.T 


2k2 


2 

3? 
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For  large  t  (when  to2 /  £2  »  1),  the  temperature  at  the  origin  is  given  by 


1  fVt  1  /  ( y2s 2 

u{0,0,0,t )  =  — 5-  /  "7t  ^ — —h 


2nlJo  {2s2  +  1)  \  £2 

l  M  l  / 

2 7T 2  -Z 0  (2s2+l) 


ds 


—s  )  ds 
£ 


1  ds+xr^-\» 


2ne  Jo  2 s2  +  1  2 TV-  2s2  +  1 


o2s2 


\/K—s^j  ds -\- T.S.T 


(14) 


i^log(2»+l)+2j«(^)+T.s.T 


where  function  q(0)  is  defined  as 


(15) 


q{%)  =  yo  2^  +  1  (M£2j2)- 


Expanding  the  term  log (2/  +  1)  on  the  righ-hand  side  of  (14),  we  arrive  at 


(16) 


aCO.Q.O.O  =  8^'°SW+  (^T«(f)  +  j^log(2))  +0(r> 


Note  that  q{a / e)  is  a  constant,  independent  of  I. 

We  now  go  back  to  the  physical  variables  tphy  and  uphy.  Recalling  (6),  we  then  find  that 


(17) 


Uphyi  0, 0,  0,  tphy)  =  «0  +  (V  0, 0,  p(K a2lphy 


For  small  values  of  tphy,  we  obtain 


Mp/ry(0,0,0,tp/,-y)  —  WoH - - - J  sft  f  1  —  -t  +  0{t  ) 

KG  27T2  V  5 


(18) 


—  Mo  H - 3  1  1 

7T  -  K~  {pCp'j  2  (J“ 
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For  large  values  of  tp/iy,  it  follows  that 


Uphyi^i  0,  tphy) 


,  20 total  &  i  r  \ 

u0+  —  i^hw) 


KG  87T£ 


KG  \2K2  V  e 


1  /  (J\  (T 

2<?  (  °  +  87T£ 


log (2)^  +  0(t  l] 


u°+  4k  K£l0g(tphy) 


(19) 


[  w 

4^R-£  ' 


2k 


4  £  /<7\  \  _ 

+  — -«(t)  )+°(^y 


{pCp)G2  J  ^KG  V  £  / 


In  the  following  we  discuss  how  to  estimate  thermal  properties  from  the  measurements  of 
the  time-dependent  temperature  at  the  center  of  the  top  surface  and  then  give  both  short  and 
long  time  behaviors  of  the  temperature  increase  at  the  center  of  the  top  surface  of  a  thin  layer. 
Statistical  errors  associated  with  parameter  estimations  will  be  evaluated  as  well. 

3.1  Estimation  of  thermal  properties  from  short  or  long  time  measurements 
of  temperature  at  the  center  of  the  top  surface 

We  consider  the  situation  where  (j)toial  (the  total  heat  flux  of  the  heat  source),  G  (the  Gaussian 
radius  of  the  heat  source)  and  £  (the  thickness  of  the  plate)  are  known,  and  upiiy(0,0,0,tphy)  is 
measured  as  a  function  of  tphy.  We  want  to  estimate  two  thermal  parameters:  pcp  (specific  heat 
capacity  per  volume)  and  K  (the  thermal  conductivity). 

The  first  step  is  to  fit  the  measured  values  of  (uphy{ 0, 0, 0,  tpi,y)  —  uq)  / ^/tphy  vs  tphy  to  a  linear 
function  of  tp\iy  for  small  tphy. 

0)  OTp/zy)  Uq 

bo  -V  tphy. 


The  coefficients  bo  and  b\  are  determined  from  measured  data  for  small  tphy  using  the  linear 
least  squares  fitting  method.  From  the  analysis  above,  we  know  that  bo  and  b\  are  related  to  the 
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unknown  parameters  as 


bo 


hotel 


3  I  , 

K1-  K2(pCj 


1  o 

)  2  CT2 


^  _  2  tytotal  tC 

3  K2K2(pcp)2(J2  ( Pcp)a 2 

Next  in  a  similar  manner  we  fit  the  measured  values  of  (uphy  (0,0,0,?^)  —  «o)  vs  log(//7/7V) 
to  a  linear  function  of  log(//7/7V)  for  large  tphy  by  the  relationship 


Uphy(fii  0;  0?  tphy)  CQ  T  C\  log(/ phy) • 


The  coefficients  c o  and  c\  are  determined  from  the  measured  data  for  large  lpi,y  in  a  linear  least 
squares  fitting  where  we  treat  log (tpiiy)  as  an  independent  variable.  Coefficients  co  and  cj  are 
related  to  the  unknown  parameters  by 

hotel 
C°  47TK-£ 


Cl 


hotel 

4  KKE 


In  turn,  the  unknown  parameters  ( pcp,  k)  can  be  estimated  from  either  (bo,b\ )  or  (co,ci). 
Suppose  coefficients  (bo,b\)  have  been  calculated  from  measured  temperature  at  short  times. 
Parameters  ( pcp ,  k)  are  expressed  in  terms  of  (bo,b\)  as 


(20) 


K  = 


3  hotel  h~irh2 
o'  3  '  uo 

z  K2G 


(21) 


Pcp  — 


2  hotel 

3  7T2  <T3 


(bob  i 


Similarly,  based  on  the  temperature  data  at  long  times,  we  can  calculate  coefficients  (co,cg 
Parameters  ( pcp,K )  can  be  written  in  terms  of  (cq,ci)  as 


(22) 


hotal  1 

4;t£  cq 


Pcp 


hotel  1  r  4  £  /a\_Cl 
2k£o2  c0exp  _y/naq\e)  cq 


(23) 
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It  is  clear  that  if  we  have  temperature  data  either  at  short  times  or  at  long  times,  we  can  estimate 
the  thermal  parameters  pcp  and  K  from  either  (20)-(21)  or  (22)-(23),  repectively. 

3.2  Number  of  unknown  parameters  that  can  be  estimated  from  short  and 
long  time  measurements  of  temperature  at  the  center  of  the  top  surface 

One  may  wonder  that  if  we  have  temperature  data  at  both  short  times  and  long  times,  would  it 
be  possible  to  estimate  more  than  2  unknown  parameters  from  the  4  coefficients  (bo ,  b\ ,  cq ,  c \ ) ? 
Below  we  show  that  it  is  not  possible  to  determine  the  parameters  if  all  three  of  ( pcp ,  K.  tytotai) 
are  unknown. 

Specifically,  we  demonstrate  that  there  are  two  different  sets  of  ( pcp ,  K.  (j)lolai  )  that  produce 
the  same  value  of  (bo,bi,co,ci).  Suppose  we  multiply  ( pcp ,  K,  <f>total )  by  a  factor  A 

(P  ^  p) new  =  A  ( P Cp ) 

Knew  —  AfC 
itytotal)  new  =  al¬ 


ii  is  straightforward  to  verify  that 


( bo)new  —  bo , 

{b\)new  —  bi, 

(pO  )new  =  6*0  ? 

(ci)new  =  C1  • 

It  becomes  apparent  that  one  can  estimate  only  two  unknown  parameters  from  the  short  time  or 
long  time  temperature  profile  at  the  center  of  the  top  surface. 

3.3  Short  time  behavior  of  the  temperature  increase  at  the  center  of  the  top 
surface 
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We  study  the  temperature  increase  (over  the  initial  temperature  liq)  at  the  origin  for  short 
time.  For  small  values  of  tpi,y,  the  leading  term  of  the  temperature  increase  takes  the  form 

uphy (0)  0, 0,  tphy)  —Uq=  — 3  —  7T  7 \frphy  U  “  7  /  VZi^phy  +  0{t  hy) 

K2K2(pcp)20I  V  3  \Pcp)° 

Thus,  the  short  time  temperature  increase  at  the  origin  satisfies  the  following  properties: 

•  It  is  proportional  to  the  total  heat  flux  of  the  heat  source. 

•  The  leading  term  is  proportional  to  the  square  root  of  time. 

•  The  leading  term  is  inversely  proportional  to  the  square  root  of  the  thermal  conductivity. 

•  The  leading  term  is  inversely  proportional  to  the  square  root  of  the  specific  heat  capacity 
per  volume. 

•  The  leading  term  is  inversely  proportional  to  the  square  of  the  Gaussian  radius  of  the 
heat  source. 

•  It  is  independent  of  the  thickness  of  the  layer. 

3.4  Long  time  behavior  of  the  temperature  increase  at  the  center  of  the  top 
surface 

Now  we  assume  a  thin  layer  and  study  the  temperature  increase  (over  the  initial  temperature 
Mo)  at  the  origin  for  long  time.  By  “thin  layer”  we  mean  e  <<  c.  That  is,  the  thickness  of 
the  layer  is  much  smaller  than  the  Gaussian  radius  of  the  heat  source.  We  first  examine  the 
asymptotic  behavior  of  function  q(£)  for  large  In  the  definition  of  q(<q)  given  in  (15),  we 
apply  a  change  of  variables  w  —  E,s  to  rewrite  q{q)  as 

1  1 

(24)  =  — —  (h(w2)  -  V^w)  dw. 

When  w  is  large,  ( h(w 2)  —  y/nw)  =  T.S.T.  Thus,  in  the  integral  above,  the  contribution  from  the 
region  of  large  w  is  negligible.  It  follows  that  for  large  q ,  in  the  region  of  dominant  contribution, 
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|t  is  small  and  we  can  expand  in  terms  of  . 


q(S)  =  r 


(25) 


y  (l— 2-^y  )  (/l(w2)  —  V/^'W;)^W;+  ■ 


<n  ,  ^3 

<3  ^ 


where  coefficients  q\  and  <73  are  given  by 


(26) 

(27) 


qi  —  {h(w2)  —  \fnw)dw, 

Jo 


q3  —  —  2w2(h(w2)  —  \/nw)dw. 

Jo 


In  Figure  3  we  plot  the  function  q(£)  and  its  approximation  ^  +  If.  As  shown  in  Figure  3, 
the  two  curves  are  getting  closer  to  each  other  when  £  becomes  bigger. 


Figure  3.  The  plot  of  #(£)  and  its  approximation. 
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Substituting  the  expansion  of  q(£,)  into  (14),  we  obtain  that  for  a  thin  layer  and  for  large 
values  of  t  the  temperature  increase  at  the  origin  is  described  by 


(28) 


tytotal 

4  KK£ 


(  log(Viy)  +  l°g  ( 


2  K 


\ 


V 


(. PCp)o 2 


4^  g4  \ 

yjn  G2  (J4  J  ' 


For  a  thin  layer,  the  long  time  temperature  increase  at  the  origin  possesses  all  of  the  following: 

•  It  is  proportional  to  the  total  heat  flux  of  the  heat  source; 

•  The  leading  term  increases  logarithmically  with  the  time; 

•  The  leading  term  is  inversely  proportional  to  the  thermal  conductivity; 

•  It  is  inversely  proportional  to  the  thickness  of  the  layer; 

•  The  leading  term  is  independent  of  the  specific  heat  capacity  per  volume  ( pcp )  and  is 
independent  of  the  Gaussian  radius  of  the  heat  source  (a);  more  specifically,  at  long 
time  tphy,  the  leading  term  of  the  temperature  increase  is  of  the  order  0(log(tpi,y))  while 
the  effect  of  pcp  and  a  on  the  temperature  increase  is  of  the  order  (9(1). 


3.5  Fitting  analytical  solution  to  measured  temperature  profile  over  a  full 
range  of  time 


We  consider  the  problem  of  fitting  the  analytical  solution  to  the  measured  temperature  over 
the  full  range  of  time.  That  is,  instead  of  fitting  to  measured  temperature  for  either  short  time 
or  for  long  time,  we  illustrate  how  to  use  all  measured  temperature  data  points  in  fitting.  As 
we  will  see  below,  by  using  the  temperature  data  over  the  full  range  of  time,  we  can  determine 
three  coefficients. 

By  combining  equations  (13)  and  (17),  we  write  the  physical  temperature  as 


(29) 


UphyiS^i  0)  tphy)  Uq 


- U 

KG 


pcpo1’ph’) 


—  Si  \f®Q}pky>  (^3); 
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where  coefficients  a\ ,  a2  and  a2  are  given  by 


cc\ 


KG 


a2 


K 

pCpO 2 


(30)  «3  =  J 

and  function  g(-)  is  defined  as 

1  fW  | 

s(w;a3)  =  ^il  W^)h{a3s]ds- 

From  measured  data  points  of  uphy(0,0,0,tpiiy)  —  uq  at  various  values  of  tppy,  we  can  estimate 
coefficients  a\,  a2  and  0C3  by  minimizing  the  distance  between  the  right  hand  side  (the  fitting 
function)  and  the  left  hand  side  (the  data).  The  theoretical  expressions  of  oti,  a2  and  a2  contain 
5  parameters: 

tytotali  ^"7  (PCp):  G 7  £• 

Once  we  know  the  values  of  «i,  OC27  and  a2,  we  can  express  k,  ( pcp )  and  e  in  terms  of  (j)totai 
and  a. 

CL\  G 


( Pcp ) 


^tytotal 
a\a2  g3 


G 

£  =  — 

0C3 

Therefore,  the  temperature  data  over  a  full  range  of  time  allows  us  to  estimate  three  parameters. 


3.6  Statistical  errors  in  parameter  estimations 


Lastly,  we  turn  to  the  problem  of  estimating  coefficients  (ai,  a2, 0C3)  from  the  measured  data 
of  upiiy(0,0,0,tphy)  at  various  values  of  tpjiy.  We  consider  the  case  of  noisy  data  where  each 
measured  value  of  uphy(0,0,0,tphy)  is 


£“^(0,0,0,^)  =  uphy(0, 0,0,tphy)  +  0.01  -N(  0, 1). 
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Here  N(0, 1)  denotes  a  random  variable  of  normal  (Gaussian)  distribution  with  mean  =  0  and 
variance  =  1.  The  measurement  noise  in  data  will  lead  to  statistical  error  in  the  estimated 
coefficients  (oq,  oq,  053).  We  study  how  the  statistical  error  varies  with  the  distribution  of  time 
instances  at  which  temperatures  are  measured.  The  result  will  serve  as  a  guide  in  selecting  time 
instances  for  measuring  temperatures  in  real  experiments. 

In  our  study  of  statistical  errors,  the  true  values  of  the  3  coefficients  are 

ai  =3,  00  =  1-5,  0(3  =  2.5. 

The  coefficients  (051,00,0(3)  are  related  to  upiiy(0,0,0,tphy)  via  function  g(w;a 3),  as  given  in 
equation  (29).  We  first  examine  the  behavior  of  function  g(w,  0(3).  Shown  in  Figure  4  are  plots 
of  g(w;  0(3)  (A)  as  a  function  of  w 2,  which  is  proportional  to  tp/jy;  (B)  as  a  function  of  w,  which 

i  1 

is  proportional  to  G/?y;  and  (C)  as  a  function  of  y/w,  which  is  proportional  to  t*h  .  It  appears 

that  the  best  way  to  robustly  capture  all  features  of  function  g(w\  0(3)  from  noisy  data  is  to  use  a 

1 

uniform  grid  in  t*h  . 

In  our  numerical  simulations,  we  use  true  values  of  (oq,  0(2, 0(3)  to  calculate  up/iy (0.0.0.! phy) 
at  various  values  of  tpiry,  and  then  add  Gaussian  noise  to  generate  a  numerical  data  set  of  mea¬ 
sured  temperatures.  Each  numerical  data  set  contains  measured  temperatures  at  1000  time 
instances.  We  compare  statistical  errors  for  3  distributions  of  the  time  instances. 

•  Distribution  1:  time  instances  uniformly  distributed  in  tpjiy  between  0  and  tfinai 

1 

•  Distribution  2:  time  instances  uniformly  distributed  in  t~hy  between  0  and  lfinaj 

1 

•  Distribution  3:  time  instances  uniformly  distributed  in  t^h  between  0  and  tfinai 
We  measure  the  statistical  error  in  (oq,  oq,  0(3)  as  follows. 


where  OClf  U'  is  the  estimated  value  of  oq  from  a  data  set  and  the  average  E  {•}  is  calculated  over 
5000  independent  data  sets.  In  Figure  5,  the  statistical  error  is  plotted  as  a  function  of  tfinai, 
respectively  for  each  of  the  3  distributions  of  the  time  instances.  When  the  time  instances  are 
uniformly  distributed  in  tphy,  the  statistical  error  is  the  largest  and  it  is  sensitive  to  the  value  of 
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(C) 

Figure  4.  Plots  of  g(w;«3)  (A)  as  a  function  of  w2,  which  is  proportional  to 

i 

tphy',  (B)  as  a  function  of  w,  which  is  proportional  to  t~h  ;  and  (C)  as  a  function 

1 

of  y/w,  which  is  proportional  to  t*h  . 

i 

t finai.  When  the  time  instances  are  uniformly  distributed  in  t*h  ,  the  statistical  error  is  the  lowest 

and  it  is  also  least  sensitive  to  the  value  of  tfmai. 

Figure  6  shows  the  statistical  error  vs  t final  for  a  large  range  of  lfmai-  It  is  clear  that  when 

1 

time  instances  are  uniformly  distributed  in  t*h  ,  the  parameter  estimation  is  the  most  robust:  it 
has  the  lowest  statistical  error  and  is  least  sensitive  to  the  value  of  t final- 

4.  Conclusions 

We  have  considered  the  temperature  distribution  for  a  semi-infinite  thin  layer  induced  by  a 
Gaussian  heat  source.  In  particular,  we  have  studied  the  asymptotic  behavior  of  the  temperature 
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Figure  5.  Statistical  error  as  a  function  of  tfinai  for  each  of  the  3  distributions 
of  the  time  instances. 


Figure  6.  Statistical  error,  for  each  of  the  3  distributions  of  the  time  instances, 
as  a  function  of  tfinat  for  a  large  range  of  tfinal. 


at  the  center  of  the  top  surface  of  the  layer.  We  have  shown  that  one  can  estimate  thermal 
properites  from  various  time  measurements  of  the  temperature  at  the  center  of  the  top  surface. 
We  have  computed  the  statistical  errors  assciated  with  the  parameter  estimations  and  found  that 
the  parameter  estimation  is  the  most  robust  if  time  instances  for  collecting  temperature  data  are 
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uniformly  distributed  in  t0  25.  As  a  future  work,  we  would  like  to  validate  our  predictions  with 
experimental  data. 
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